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Abstract. We investigate the low-temperature reaction rates for radiative capture processes of three par- 
ticles. We compare direct and sequential capture mechanisms and rates using realistic phenomenological 
parametrizations of the corresponding photodissociation cross sections. Energy conservation prohibits se- 
quential capture for energies smaller than that of the intermediate two-body structure. A finite width 
or a finite temperature allows this capture mechanism. We study generic effects of positions and widths 
of two- and three-body resonances for very low temperatures. We focus on nuclear reactions relevant for 
astrophysics, and we illustrate with realistic estimates for the a-a-a and a-a-n radiative capture processes. 
The direct capture mechanism leads to reaction rates which for temperatures smaller than 0.1 GK can be 
several orders of magnitude larger than those of the NACRE compilation. 



PACS. 21.45.-v Few-body systems - 25.40.Lw Radiative capture - 26.20.-f Hydrostatic stellar nucleosyn- 
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1 Introduction 

The techniques employed to solve few-body problems have 
improved substantially over the last decade. This is not 
only due to the computer development but also to the 
higher efficiency of both the theoretical formulations and 
1 the numerical methods. In particular, a number of three- 
body problems are now all solvable at least within a cer- 
tain requirement for the accuracy. The easiest of these 
problems are bound states [312] and perhaps resonances 
[3,4,5,6,7,8 , whereas the more challenging are continuum 
■ structures in general p fTUlITT] . The most demanding are 
| reactions where the structures before and after the process 
are different [12l[T3lfT4] . 

Applications to nuclear astrophysics have been in de- 
mand for years [TMTHlfT7t fTS l [T ^ [^D ] . In particular, radia- 
tive nuclear three-body capture processes are indispens- 
able in bridging the A = 5,8 gaps [TTIl2"Tll2"2II23| . but also 
building up the heavier nuclei may involve three-body re- 
actions. Also, to pass the waiting points in the rapid pro- 
ton capture process necessarily involves three particles, 
i.e. two protons and a core [Mll2"5"llT§] . Along or outside 
the neutron dripline, or perhaps along the rapid neutron 
capture process, three particles may be involved as well in 
decay, capture, or reaction processes [T2, 8,26,2711128]. 

Traditionally, that is mostly for convenience of com- 
putation, the three-body reactions are described as two 
subsequent two-body reactions ^5lf2TlfT6ir29] fT9 ] [30 ] . This 
two-step sequential mechanism has proved to be rather ef- 
ficient in most cases but it cannot always be accurate. For 



example, an intermediate configuration of some stability 
has to be present, and resonances offer themselves as these 
stepping stones. However, they may not be present or they 
may be too unstable, such that they decay too quickly af- 
ter they are formed. It is also well-known that important 
contributions arise from the continuum background to the 
resonances j5| l201lllj . This is often called the direct part. 
All such difficulties with the sequential models disappear 
by using the full three-body formulation of the reactions. 
Then the sequential mechanisms are still possible but now 
as a contributing process which does not have to be treated 
separately [3T1I5U] . 

In reactions in stellar environments the temperature, 
rather than the energy, is often the parameter. For high 
temperature, if one of the two-body subsystems shows a 
relatively low-lying narrow resonance, the process is fairly 
accurately described by sequential model. In this case, al- 
though the three-body background contribution is miss- 
ing in the description of the process, it is either small or 
appropriate parameter adjustments account for this part. 
When the intermediate two-body state is missing (or it is 
too broad), the strength function is continuous, without 
a significant peak structure arising from resonances |121 
|2"8] , and the reactions can be described fairly well with 
knowledge of the three-body continuum and appropriate 
energy averaging [12„20"1[2"5]. 

At small temperature, where important astrophysical 
reactions take place, only very low energies are involved, 
which are typically clearly smaller than the energy of the 
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possible intermediate two-body states. As a consequence, 
the situation concerning the reaction mechanism is un- 
clear. In any case, the complete computation of the di- 
rect process in coordinate space is very inaccurate due to 
the fact that a narrow grid of discrete continuum states 
is required in the very low energy region. This is often 
unaffordable from the numerical point of view, and the 
problem is usually achieved within the formalism by as- 
suming a specific form of the strength function or equiv- 
alently of the photodissociation cross section 12 , 20 , 28 . 
This has the advantage of being a measurable quantity, 
and thus open to tests. To circumvent through pure theo- 
retical calculations would require an accurate model and 
many points at small energies above the three-body thresh- 
old. 

In view of the improved three-body techniques it is 
surprising that the sequential model still has not been 
systematically tested against genuine three-body compu- 
tations. Even though such three-body results for small en- 
ergies above threshold are inaccurate, it is still possible to 
compare direct and sequential rates by using the different 
model formulations with the same common input. This 
is the purpose of the present paper which focuses on the 
low temperature region such that only three-body ener- 
gies below the intermediate two-body resonance energies 
are relevant for the reaction rate. For these temperatures 
reliability of the sequential model is more than dubious. 
The idea is to take the established and presently applied 
formulations for sequential capture, and compare to a sim- 
ilar formulation which includes the three-body direct cap- 
ture. These calculations are based on the assumption of 
a specific form of the strength function (or the photodis- 
sociation cross section), and therefore do not require the 
precise numerical solution of the complicated three-body 
problem. Even though this is a drastic simplification of 
the problem, a comparison of the results obtained with 
the two possible mechanisms will permit to enlighten the 
importance of the direct capture mechanism at low tem- 
perature. In section II we give the appropriate definitions 
and formulae for later use. In section III we compute rates 
in the direct and sequential models and compare results 
as a function of the position and the widths of the two and 
three-body resonances. In section IV and V we apply the 
general method to the triple alpha process and the similar 
radiative capture leading to 9 Be. Finally we summarize 
and conclude in section VI. 



2 Direct versus sequential rates 

A method computing the reaction rates for a three-body 
radiative capture reaction without assumptions about the 
underlying mechanism was employed previously in |12,20, 
|2"5] . In the first part of this section we give a simplified de- 
scription of this method, which, in particular, should be 
used when a direct capture takes place. However, when de- 
scribing these processes, the corresponding reaction rates 
are traditionally described by a sequential mechanism via 
an intermediate structure [To1[T51[2^1[2Tj . The main aspects 
of this description are summarized in the second part of 



this section. By use of the different expressions given, a nu- 
merical comparison between both descriptions is allowed, 
in particular on the very low temperature region. This is 
done by use of the energy dependent penetration factors 
given in the last part of this section for each of the two 
possible capture mechanisms. 



2.1 The direct picture 

A three-body radiative capture process is usually denoted 
by a + & + c — > A + j, and the reversed reaction is A + 
7 — > a + b + c, where a, b and c denote the three particles. 
The reaction rate at a given three-body kinetic energy E, 
Rabc(E), can be obtained by use of Eqs.(20) and (32) in 
Ref. [T5], which lead to: 



Rabc(E) 



h 3 



2g, 



C 2 (HabfJ-ab,c) 3/2 3a9b9c 



(i) 



where E 1 = E + \B\ is the photon energy, B (< 0) is the 
three-body energy of the nucleus -A, ix a b ctnd fj>ab,c &re the 
reduced masses of the a-b two-body system (related to the 
Jacobi coordinate x) and the ab-c system (related to the 
Jacobi coordinate y), respectively [2J, gi = 2J* + 1, where 
Ji is the angular momentum of particle i (i — a, b, c, A), 
and v is the number of identical particles in the three-body 
system. Finally, a 1 {E 1 ) is the photodissociation cross sec- 
tion of the nucleus A. 

The energy averaged reaction rate is obtained as a 
function of the temperature T by using the Maxwell-Boltzmann 
distribution as weighting function. For three-particles the 
Maxwell-Boltzmann distribution takes the form: 



1 E 2 



(2) 



where ks is the Boltzmann constant. We then obtain the 
following expression for the energy averaged reaction rate 



(Rabc(E)) = v\ 



fi 3 



8tt 



9 A 



c 2 {nabLiab, c ) 3/2 g a gbg c 
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E^a^E^e "B^dE 



(3) 



Therefore, once the photodissociation cross section cr 7 for 
the process A + -f — >• a + b + c is known, the rate (R a bc{E)) 
can be easily obtained. For instance, it can be obtained 
directly from the experimental <7 7 cross section. 

Numerical calculation of a 1 for the A + j — > a + b + c 
process is a delicate task, especially at low temperatures. 
For instance, the method used in [T2"ll2"01l2"8] would require 
an extremely narrow grid of discrete continuum states. 
To achieve this, the continuum wave functions should be 
computed in a box of a sufficiently big size, which in turn 
would imply the need of an enormous basis set. All in 
all, the method becomes inefficient for such low temper- 
atures. However, when the photodissociation proceeds by 
populating a three-body Breit-Wigner shaped resonance 
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of particles a, b, and c with total angular momentum J 
and energy En, then a^(E^) can be written as: 



£7 7 (£ 7 ) 



2 J 



r abc (E)r 7 (E) 



2 9A ei (E-E R y + r{Ey/e 



(4) 



where T 7 and r abc are the partial decay widths of the reso- 
nance for gamma and three particle emission respectively, 
and r — r a b c + r y is the total width. This is analogous 
to the method used in the usual sequential picture that is 
described in the next section. In this way, by using Eq. © 
the results provided by the two different schemes can be 
easily compared. 

It is important to note that Eq.© has been derived 
without making any assumption concerning the reaction 
mechanism leading to the formation of the nucleus A. It is 
therefore completely general. In particular, this is the ex- 
pression to be used when the capture mechanism is direct, 
meaning that the capture takes place without populating 
any intermediate two-body state. In this case the width 
r a bc{E) corresponds to the width for direct decay of the 
three-body resonance into particles a, b, and c. For this 
reason we shall refer to Eq.© as the reaction rate in the 
"direct" or "three-body" picture. 



2.2 The sequential picture 

Let us now assume that one of the internal two-body sub- 
systems, for instance the a-b system, shows a relatively 
narrow two-body resonance at some energy E r . The usual 
procedure to understand reactions of the type a + b + c — > 
A + 7 is then to interpret it as two sequential two-body 
processes [11)1121)] . In the first step particle a captures b to 
populate the intermediate a-b two-body resonant state. In 
the second step, the a-b system is able, before decaying, 
to capture particle c, populate some three-body resonance 
of the nucleus A, and then decay by photo emission into 
one of the bound states of A. 

The reaction rate for such a two-step process is given 
by the rate for the capture of c by the two-body subsystem 
a-b, (R a b,c{E" , E')), weighted with the rate for formation 
of a-b 161: 



(R abc {E",E')) = 



A 8irh 



(J-ab 



3/2 



1 + S ab fi- ab \2nk B T / 
r ^^e- E ''/ k * T {R a bAE",E'))E"dE", 

JO l ab\& ) 



(5) 



where the total three-body energy, E = E' + E" , is given 
in terms of the relative energy, E" , between particles a and 
b and the energy, E' , of particle c relative to the center of 
mass of a-b. The function 8 ab is 1 if a and b are identical 
particles, and otherwise. 

The elastic a-b cross section, a a b, in the equation above 
takes the form: 



<7ab{E") = (1 + S ab 



Jab 7T 



r ab (E" 



g a g b ^ (E" - E r y + {r ab {E")Y/A 



, (6) 



where k 2 = 1^ a bE" /ti 2 , E r is the energy of the resonance 
in the a-b system and r a b(E") is the corresponding width, 
and g a b = 2J a b + l, with J a b being the angular momentum 
of the two-body resonance. 

Also, following Ref.[16], we have that 



(Rab,c{E" , E')) — — 

^ab,c 



^ab,c 



3/2 



2 7 rfc B T / 

POO 

x / a ab AE",E')e- E 'l k - T E'dE', (7) 
Jo 

where the cross section cr a b.c{E" ,E') for the capture of 
particle c by the two-body subsystem a-b is related to 
the photodissociation cross section through the detailed- 
balance (or reciprocity) theorem, which leads to |29j : 



where E~ 



<TabA E ') 



E+ \B\ 



9 A 



, £°~f(E 7 ) , (8) 

g a b9c fJ-ab,c C Z E' 



Similarly to the direct process, the photodissociation 
cross section tr 7 again takes the form: 

. 2 J + 1 vr^c 2 r ab , c (E')r 7 (E' + E") 



2g, 



E 2 (E-E R ) 2 + r(E',E") 2 /A- 



This expression is formally identical to Eq.@, but now 
Eab,c refers explicitly to the partial width for decay of the 
three-body resonance with angular momentum J into the 
two-body resonance a-b plus particle c. The assumption 
is that no other decay mode exists, and the direct decay 
circumventing this two-body path is not allowed, or at 
least negligibly small. As before, _T 7 is then the partial 
width for gamma decay, and r = r a b >c + J" 7 is the total 
width. 

Replacement of Eq.© into Q, and of Eq.© into © 
leads then to the following expression for the energy av- 
eraged reaction rate: 



(R abc (E",E')) 



<JA 



SPi 



1 



5ab 9ab9c irc 2 (k B T) 3 n 1 ^^ 1 , 2 



E"^^dE» T E^(E,)e-WdE (10) 

1 ab{& ) JE" 



where we have replaced the dependence on E' by E by 
use of E = E' + E". We shall refer to the expression in 
Eq.(frU]) as the reaction rate in the "sequential" picture. 



2.2.1 Extreme sequential limits 

In this subsection we shall derive the form that the reac- 
tion rate (|T0l) takes in the extreme cases of an infinitely 
narrow three-body resonance and/or an infinitely narrow 
intermediate two-body state. We shall refer to these rates 
as the ones obtained in the "extreme sequential" picture. 

The connection between the reaction rates © and (|TU| 
is made evident for the particular case of a very narrow 
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a-b two-body resonance (r ab very small compared to the 
energy of the resonance). In this case we have that Eq.© 
can be replaced by: 

Vgb(E") , g ab 2ir 2 „ 

from which Eq.([TU]) becomes: 



(R abc (E)) = v\ 



h 3 8tt 



9A 



C 2 (fJ.ab^ab,c) 3/2 9a9b9c 



E 2 a 7 (E 7 )e k B^dE. (12) 



(k B TfJ Er 

and the three-body energy E is given by E r + E' . 

We can immediately see that the reaction rates in 
Eqs.© and (TT21 are identical except for the lower limit 
in the integration. From in the first case and from the 
intermediate two-body resonance energy E r in the sec- 
ond. This reflects the fact that in the extreme sequential 
picture described by Eq. (|12[) . where an infinitely narrow 
intermediate two-body resonance is assumed to be popu- 
lated, the total three-body energy E must be larger than 
E r . For three-body energies smaller than E r this sequen- 
tial picture cannot provide any decay, since the energy 
is too small to populate the two-body intermediate reso- 
nance. This means that in this extreme limit, when the 
three-body resonance energy En is smaller than E r , the 
sequential decay has to proceed through the tail of the 
cross section in Eq.© corresponding to E > E r > Er. 

We can also consider the case where the width r abc 
in Eq.© is very small. Assuming that _T 7 <gc r abtC is still 
fulfilled, we then have from Eq.© that: 



Vj(Ej) 



which leads to 



2j + ir^^ 



fjA 



g a g b gc (k B ± r {^ab^ab,c) A/z 



E n 



r ab (E") 



(E» - E r ) 2 + {r ab {E")) 2 /A 



dE". 



(14) 



From the expression above we can see that when En < E r 
the decay in this limit must proceed through the tail of the 
cross section in Eq.© corresponding to E" < En < E r . 

When the widths -T a h,c and r ab are both very small 
then the approximations (jlip and p3[) can be used simul- 
taneously, and either Eg . (fT2"|) or (TT4")) lead to the following 
expression for the energy averaged reaction rate: 

g a g b gc (kb^t [Hab^abA 1 

which is valid when the two-body resonance energy E r is 
below the three-body resonance energy En- Otherwise, if 
En < E r the reaction rate is zero in this limit of very 
small widths. The expression in Eq.(fT5]) does not depend 



on the energy E r of the intermediate two-body resonance, 
and it is actually the same expression that one gets when 
the approximation (|13[) is used on the reaction rate Q 
obtained in the direct picture. 

At this point it is important to remind ourselves that 
Eqs. (fl~2"l) . (JH) , an d dSl) are limiting cases involving very, 
or even infinitely, narrow resonances. In particular, Eg. (1121) 
permits to see very easily the connection between Eqs.@ 
and (fTUjl . However, correct sequential calculations have to 
be performed by use of Eo. lfTO)) . as done for instance in 
PH or [IS]. 



2.3 Low energy dependence 

The Boltzmann exponents in Eqs.Q and (p~0|) only allow 
contribution from three-body energies up to a few times 
ksT. For instance, for a typical temperature in the core 
of a star of 10 7 K we have that k B T « 1CT 3 MeV, and 
only energies of the order of a few keV's and below can 
contribute to the rate. These energies are very small at the 
nuclear scale, where excitation energies are typically of the 
order of the MeV. Therefore, for such low temperatures, 
the reaction rate in Eq.Q or (fTO)) is dictated by the low 
energy tail of the corresponding cross sections. 

As seen in Eqs.flU, ©, and © the behavior of the 
cross sections at very low energies is determined by the 
energy dependence of the partial widths, which are pro- 
portional to the penetration factor through the barrier re- 
sponsible for the corresponding resonance (r(E) oc P(E)). 
The constant of proportionality is determined by assum- 
ing that when evaluated at the resonance energy 
width r(E res ) provides the established (or experimental) 
width of the resonance 2"b . This is therefore leading to the 
expression: 

where the explicit form of the penetration factor P(E) 
should be determined for each particular case. 

For gamma decay it is well known from text books (see 
for instance the appendix E.l of [35]) that the decay con- 
stant, and therefore the penetration factor and r-y(E), is 
proportional to E 2X+1 : where E 7 = E + \B\ is the pho- 
ton energy and A is the multipolarity of the electromag- 
netic transition. As discussed in Eq. (ll6[) . the proportion- 
ality constant is determined in such a way that r i {En) is 
equal to the experimental value -T 7 for the gamma decay 
width. The parametrization is then: 



r 1 (E) = r, 



( E+\B\ 
\E R + \B\ 



2A+1 



(17) 



In the same way, when the Coulomb interaction is 
not involved, the partial widths r abiC (E') and r ab (E") in 
Eqs.© and (0) for the sequential picture take the form: 



E a b,c(E') — r a (, 



E' 



En — E r 



£ ab .a + l/2 



E' 

r ab {E") = r ab ( -g- 



(18) 
(19) 
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where £ a b,c and l a b are the relative orbital angular mo- 
menta between particle c and the center of mass of a-b 
and between particles a and b, respectively, and where 
we have taken into account that the energy dependence 
of these partial widths is given by Eqs.(22) and (23) of 
Ref.[33]. 

When the Coulomb interaction is decisive, insertion of 
Eqs.(lO) and (11) into (3) of Ref.[33] leads to the following 
energy dependence of the partial widths in the sequential 
picture: 



r a b,c{E') — r a b. c 
r a b{E") = r a b 



f e 2b ab ^/^E T1 -E r 
1 _|_ e 1h ab ^l-JW 
1 + e 2WV^ 



1 + e 2b ab /VE*> 



where 



b a b,c = -^(Z a + Z b )Z c e/ 



2fiab,< 



ZnZhe 1 



%Hab 



(20) 
(21) 

(22) 
(23) 



and where Z a , Zb, and Z c are the charges of particles a, 
b, and c, respectively, and e is the electron charge. 

For direct decay we get, from Eqs.(12) and (3) of [33], 
for Coulomb potentials that: 



r a bc{E) — r a 



1 



be " 



1 -I- e 2b abc /VE 



(24) 



with 



babe = 



h 2 (m a + mb + m c ) 



1 ' ZiZje 2 ) 2 ^ 3 (niim.j 



3/2 



(25) 

where m a , nib, and m c are the masses of the particles and 
the sum runs over the three possible pairs of particles. 

For direct decay, without Coulomb potentials, we get 
the energy dependence of r a b c (E) in Eqs.Q and ^ from 
Eq.(21) of Ref.[33]: 



r a bc(E) — r a 



be 



E 
~E~n 



: +£ ab +2 



(26) 



3 Parameter dependence 

The computed reaction rates depend on the choice for 
the reaction mechanism, i.e. either direct or sequential. 
To illustrate it we start with a system of three identical 
charged particles. We shall assume that this three-body 
system has a resonance at the energy Eji — 0.15 MeV 
and -T 7 = 10~ 9 MeV. All the degeneracy factors gi are 
taken equal to 1, the charges are taken twice the proton 
charge, and the masses equal to the mass of the alpha 
particle. 

If the capture process is assumed to proceed directly, 
the reaction rate is then given by Eqs.Q and ([3]), and 



for three charged particles, the energy dependence of the 
width for particle decay is given by Eas. ([24"]) and (j25|) . 
On the other hand, if we assume that the capture process 
takes place sequentially through an intermediate two-body 
resonance, the reaction rate is instead given by Eqs. ([5]), 
©, and (ITU1) . and the energy dependence of the widths, 
r a b and r a b tC , for particle decay in Eqs.© and ([9]) are 
now as in and ([2TJ)) . 

As a first case, let us assume that the three-body and 
the intermediate two-body resonances are both narrow, 
e.g. by choosing r fl = 10~ 5 MeV and T r = 10~ 6 MeV. 
The computed reaction rate in the direct picture is then 
given by the thick solid line in the upper part of Fig.[T] For 
the sequential case, when the capture proceeds through 
the intermediate two-body resonance, as seen in Eqs.© 
and (fTTJ|) . the reaction rate depends as well on the energy 
of the two-body resonance. We have chosen the energies 
of 0.05 MeV, 0.10 MeV, 0.13 MeV, and 0.20 MeV, as indi- 
cated by the straight lines in the right part of Fig [2j The 
straight line in the left part indicates the energy chosen 
for the three-body resonance. In this case of very nar- 
row three-body and two-body resonances, the lorentzians 
describing them can not be distinguished from the lines 
indicating the energies in Fig. [5] The corresponding reac- 
tion rates in the sequential picture are given by the short- 
dashed, dot-dashed, and long-dashed curves in the upper 
part of Fig. [T] for the three cases where E r takes the values 
0.05 MeV, 0.10 MeV, and 0.13 MeV, respectively. 

As we can see in the upper part of Fig. [TJ for tem- 
peratures below about 0.02 GK the reaction rates in the 
two pictures differ from each other. In the first two cases, 
E r = 0.05 MeV and E r = 0.10 MeV, the sequential rate 
is below the direct one, while for E r — 0.13 MeV, the se- 
quential rate at low temperatures is above the direct one. 
This behavior is determined by the interplay between r a b jC 
and r a b in Eas. ([2"0|) and (|2"Tj) . When E r increases and ap- 
proaches Er, then r a b. c increases, implying a more signif- 
icant tail in the cross section ©, and therefore favoring 
an increase of the reaction rate at low temperatures. How- 
ever, at the same time the increase of E r decreases r a b, 
which for the same reason tends to reduce the reaction 
rate. The final result of the competition between these 
two effects pushing in opposite directions depends on the 
masses and charges of the particles involved in the pro- 
cess. In our case the total effect amounts to an increase 
of the reaction rate when E r increases. The comparison 
between the direct and the sequential calculations is bet- 
ter seen in the inner part of the figure, where the different 
curves show the ratio between the rates in the direct and 
sequential pictures. As we can see, the difference between 
them reaches up to 12 orders of magnitude for T = 0.001 
GK. 

For temperatures higher than about 0.03 GK the rates 
obtained in the sequential picture are not sensitive to the 
properties of the intermediate two-body resonance, and 
the computed rate agrees with the one obtained in the di- 
rect picture. This is consistent with the fact that when the 
three-body resonance and the two-body one are both very 
narrow the reaction rate can be approximated by Ea. (|15p . 
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Fig. 1. (color online) Reaction rates for the radiative capture 
of three identical charged particles (with charge and mass equal 
to the ones of the a-particle). The widths Er and E r of the 
three- and two-body resonances are 10 -5 MeV and 10 -6 MeV, 
respectively in the upper part, 10 -5 MeV and 0.1 MeV, respec- 
tively in the central part, and 0.1 MeV and 10~ 6 MeV, respec- 
tively in the lower part. The thick-solid curve shows the reac- 
tion rate in the direct picture. The short-dashed, dot-dashed, 
and long-dashed curves are the rates in the sequential picture 
for two-body energies E r = 0.05 MeV, 0.10 MeV, and 0.13 
MeV. The ratio between the rates in the direct and sequential 
pictures are shown in the insets. The thin solid line is the rate 
in the extreme sequential limit. The dotted curve is the reac- 
tion rate obtained when E r = 0.20 MeV (Er < E r ) (see text). 
Na is the Avogadro's number. 
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Fig. 2. (color online) Energy diagram showing the three-body 
resonance energy Er and two-body resonance energies E r used 
in the schematic calculations in section III. The curves are the 
lorentzians corresponding the resonances when having a width 
of 0.1 MeV. 




I I I I I I | 1 1 — I I I I I 





=0.15 MeV, r R = 


10 


5 MeV 


- E r= 


=0.15 MeV, r R = 


10 


15 MeV 


..... Er= 


=0.15 MeV, r R = 


10 


25 MeV 


E R= 


=0.15 MeV, r R = 








10" 
T(GK) 



10 



10 



Fig. 3. (color online) Reaction rates in the sequential picture 
for the same system as in the upper part of Fig. [T]with E r = 
0.10 MeV but for different small values for the width of the 
three-body resonance Er. The thick-solid and thin-solid curves 
are the same as the dot-dashed and thin-solid curves in the 
upper part of Fig. [T] The dashed and dot-dashed curves are 
calculations with Er = 10" 15 MeV and Er = 10" 25 MeV, 
respectively. 



which depends neither on E r nor on r r . This is referred 
to as the rate in the extreme sequential picture, and it 
is shown in the upper part of Fig. [T] by the thin solid 
line. As we can see, for T > 0.03 GK this approximation 
agrees with both the full direct and the sequential cal- 
culations. At low temperatures this is not true anymore, 
which shows that the finite width of the resonances in- 
fluences the rates, and the limit of zero width leading to 
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Eg. (|15[) becomes increasingly invalid. The reason is that 
the effect from the tails of the resonances, where the pene- 
tration factors play a role, is completely absent in Eq. (fl~5")) . 
This is easily recognized in Fig. [21 where the thick-solid 
and thin-solid curves are the same as the dot-dashed and 
thin-solid ones in the upper part of Fig. Q] (sequential reac- 
tion rate with E r — 0.10 MeV and reaction rate in the ex- 
treme sequential limit, respectively), while the dashed and 
dot-dashed curves are the calculations with Tr = 10~ 15 
MeV and = 10~ 25 MeV, respectively. As we can see, 
the narrower the three-body resonance the lower the tem- 
perature at which the full sequential calculation matches 
the curve in the extreme sequential picture. 

The central and lower parts of Fig. [1] show the same 
reaction rates as in the upper part, but for r r — 0.1 MeV 
and r*R = 0.1 MeV, respectively. A feeling of the overlap 
between the three-body and two-body resonances in the 
different cases can be got from Fig. [5J where the curves 
plotted are lorentzians with r = 0.1 MeV. The behavior 
observed in the central and lower parts of Fig. [1] is similar 
to the one found in the upper part, where both resonances 
are simultaneously narrow. Also the ratio between the di- 
rect and sequential calculations (inner part of the figures) 
is similar. In these two cases, the extreme sequential curves 
correspond to the ones obtained with the expressions (|14l) 
and (dU), respectively. As we can see, Eq. p^tj) does actu- 
ally depend on the properties of the two-body resonance 
energy, which leads in this case (central part of Fig. [TJ) to 
some small differences between the different calculations. 
In particular, the curve in the figure (thin solid line) has 
been obtained using E r = 0.10 MeV. 

It is important to note that the equations in Ref . |33j 
which have led to Eqs. flTS)) and (f2~Tj) have been obtained 
using the adiabatic approximation to describe the three- 
body system. The penetration factors are then extracted 
from the WKB tunneling probability through the poten- 
tial barrier of the lowest adiabatic potential. When one of 
the two-body systems has a resonance, this lowest adia- 
batic potential shows a flat region at the two-body reso- 
nance energy, and eventually goes to zero. The narrower 
the two-body resonance, the longer its lifetime, and the 
larger the flat region in the adiabatic potential. In fact, 
in the limit of zero width (or infinite lifetime) we actually 
have a bound two-body state, and the adiabatic potential 
goes asymptotically to the two-body energy of the bound 
state. This implies that the expressions (|T5|) and (j2T))) are 
meaningful only when the three-body resonance energy 
Er is bigger than the two-body resonance energy E r . 

In the opposite case, when Er < E r , the reaction 
mechanism is dominantly direct, since the intermediate 
configuration is energetically forbidden. However, it is still 
possible to exploit the sequential mechanism but obviously 
not until the decay is completed. The distance between 
particles must be limited, or equivalently the time allowed 
for the non-energy conserving process is limited. If the de- 
cay still has to take place through the sequential decay 
mechanism then the two-body resonance width has to be 
finite (see the dot-dashed curve in Fig. [5J where even if 
it corresponds to Er < E r a significant part of the reso- 



nance tail falls in the energetically allowed region). This 
mechanism was previously named virtual sequential decay 

mi- 

From the discussion above we then have that, when 
Er < E r , the barrier determining the penetration factor 
is much bigger than in the opposite case, and in fact the 
narrower the two-body resonance the thicker the barrier. 
As a consequence, for Er < E r we have that Eqs. ([!%)) and 
d20J) can not be used, and the three-body width is expected 
to be much smaller than when Er > E r . For this reason, 
the reaction rate for the cases where Er < E r can only be 
estimated from Eq. (fT4]) . where the value of r a f, tC has been 
assumed to be zero. From this equation, since Er < E r , 
we can immediately see that the reaction rate is produced 
by the tail of the intermediate two-body resonance (dot- 
dashed curve in Fig. [5]). Clearly, the broader this resonance 
the larger the reaction rate. This is shown by the dotted 
curves in the upper and central parts of Fig. [T] for r r = 
10~ 6 MeV and T r = 0.1 MeV, respectively. These two 
calculations have been made with resonance energies E r = 
0.20 MeV and E R = 0.15 MeV. 

In the central part the width of the two-body resonance 
is not small, E r =0.1 MeV, and the difference with the 
previous calculations is hardly seen in logarithmic scale, 
although it is two orders of magnitude smaller than the 
direct picture calculation. When F r = 10~ 6 (upper part), 
the dotted curve is up to eight orders of magnitude smaller 
than the other calculations, consistent with the fact that 
in the limit where r r = the rate provided by Eq. (fT2"|) is 
zero. In the lower part of the figure the corresponding cal- 
culation with E r = 0.20 MeV is not shown. First, because 
for Er < E r the width of the three-body resonance is ex- 
pected to be very small, which is not the case in this cal- 
culation, and second, because even if we insist in comput- 
ing it, the penetration factors (fT5)) and (I2U1) are obviously 
meaningless for Er < E r , and another parametrization 
would have to be designed. 

For completeness, we show in Figs. [3] and [3] the same 
calculations as in Fig. [1] but only for the cases where, 
respectively, two or none of the particles forming the two- 
body resonance are charged. The meaning of the curves 
is as in Fig. [TJ and the only difference lies in the use 
of the proper penetration factors for the cases where the 
Coulomb interaction can be neglected. 

In particular, in Fig. 2] we assume that the two charged 
particles are the ones forming the intermediate two-body 
resonance in the sequential picture. Therefore, the pene- 
tration factors to be used are ([M]) in the direct case, and 
(|T8| and (j2Tj) in the sequential case. In this case, when E r 
increases and approaches Er we have that -T a b iC increases 
as a polynomial, while r a b decreases exponentially. As a 
consequence, the decrease in the reaction rate due to the 
decrease of r a b dominates over the increase produced by 
the increase of -T &,c, and the final result is that the com- 
puted reaction rate decreases at low temperatures when 
E r increases, although the difference between the calcula- 
tions when E r — 0.10 MeV and E r = 0.13 MeV is hardly 
seen in the figure. In Fig. [5] the penetration factors to be 
used are (121))) in the direct case, and (fT5]l and (|T5)) in the 
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Fig. 4. (color online) The same as in Fig. [1] when only the 
two particles forming the intermediate two-body resonance are 
charged. The centrifugal barriers correspond to orbital angular 
momentum £ a b,c = 2. 
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Fig. 5. (color online) The same as in Fig. [T]for three identi- 
cal uncharged particles. The centrifugal barriers correspond to 
orbital angular momentum £ a t,c ~ 2. 



sequential case. The general trends of the reaction rates 
shown in Figs. [4] and [5] are similar to the ones shown and 
discussed in Fig. Q] The main difference is the absolute 
value of the rates at small temperatures, which increases 
enormously when the Coulomb barrier disappears. 

These general features of the reaction rates, especially 
the ones concerning their dependence on model assump- 
tions of the capture process, can now be used in practical 
cases of (astro)physical interest. We shall in the next two 
sections focus on the triple alpha reaction (three identical 
charged particles with narrow three-body and intermedi- 
ate two-body resonances), and on the a + a + n — > 9 Be + 7 



process (two identical charged particles forming a narrow 
two-body resonance, and a relatively broad three-body 
resonance). 



4 The triple alpha rate 

The reaction rate for the triple alpha reaction, a+a+a — > 
12 C + 7, at very low temperatures is quite controversial. In 
fact, in the recent work [33], the reaction rate at a temper- 
ature of 10 -2 GK has been found to be about 20 orders 
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of magnitude bigger than the one given in the NACRE 
compilation |16j . Such enormous increase in the reaction 
rate would have dramatic consequences for the late stages 
of the stellar evolution in low mass stars [35] . 

As demonstrated in the previous sections, the key quan- 
tity to obtain an accurate reaction rate for a radiative cap- 
ture process is the photodissociation cross section for the 
inverse reaction. This means that in our case a 1 for the 
12 C+7 — > a+a+a reaction is needed. This cross section is, 
to a large extent, determined by the resonance spectrum 
in 12 C, which is rather well known except for the lowest 
2 + resonance, whose excitation energy is not well estab- 
lished yet. Fortunately, as shown in [50], this resonance is 
playing a role in the reaction rate only for temperatures 
higher than about 2 GK. Therefore, this uncertainty in 
the 12 C spectrum is irrelevant for the analysis of the re- 
action rate at very low temperatures, which is dominated 
by the E2 transition + — s> 2^ 2(3 . This notation defines 
the transitions from the + continuum states in the three 
a system to the bound 2 + state in 12 C. 

Furthermore, as also shown in [20] . for temperatures 
higher than 0.1 GK the computed contribution to the re- 
action rate from the + — > 2^ transitions is very similar 
when performing a full three-body calculation and when 
assuming a sequential capture through the low-lying + 
state in 8 Be (0.092 MeV above the two-body threshold). 
This is due to the existence of the narrow + Hoyle three- 
body resonance in 12 C at about 0.38 MeV above thresh- 
old. This state heavily dominates the calculation, which 
has a large strength corresponding to strong population 
and subsequent decay of the Hoyle state. Three-body cal- 
culations have shown that the particle decay of the Hoyle 
state proceeds almost fully sequentially [35] , and therefore 
so does also the inverse process. 

However, for very small temperatures the relevant three- 
body energies in the reaction rate are clearly below the 
0.092 MeV of the intermediate + resonance in 8 Be. It is 
then not so obvious that the sequential picture through 
that + state is still appropriate. The occurrence of a di- 
rect or a sequential capture mechanism would imply a dif- 
ferent behavior of the tail of the photodissociation cross 
section (see Eqs. ([2"0"|) and ([2"3])). and therefore it would lead 
to a different reaction rate. This is in fact observed in the 
schematic case in the upper part of Fig. [1] where an ex- 
ample similar to the triple alpha reaction showed a model 
dependent rate varying by several orders of magnitude at 
low temperatures. 



4.1 The photodissociation cross section 

As mentioned above, for temperatures higher than 0.1 GK 
the direct and sequential descriptions provide similar re- 
sults. Therefore, for T > 0.1 GK, the simplest choice for 
the photodissociation cross section of 12 C is the sequential 
expression given by Eq. © . We use the same values for the 
energies and widths of the different states in 12 C and 8 Be 
as the ones specified in [TB] . The energy dependence of the 
widths is then as given by Eq. ([T7jl (with A = 2) for -T 7 (£), 
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Fig. 6. (color online) Photodissociation cross section as given 
in Eq.([9]l for the 12 C+7 — > a+a+a reaction as a function of the 
three-body energy E. The parameters are as given in [X6| . The 
solid line has been obtained assuming a sequential mechanism, 
using the widths as given in (|20|l and (|22p . The dashed line 
corresponds to a direct mechanism, using the widths as given 
in (|24|) and (|25|) . The inner part is a zoom of the cross section 
in the matching region. 

and as given by Eq. (l2Tfl) (with the 6-parameter (|2"2"10 for 

r a b,c(E). 

If we now assume that the same energy dependence, 
corresponding to a sequential mechanism, is still valid for 
low energies, we then get the photodissociation cross sec- 
tion shown by the solid line in Fig. [6j The peak at about 
0.38 MeV corresponds to the Hoyle state. 

The other possibility is that at low energies the pho- 
todissociation cross section be dominated by a direct pro- 
cess. In this case the tail of the cross section would be given 
by Eo. flMl) . where the 6-parameter is now as in Eq. (12~5|) . 
Then Eq.Q leads to the cross section shown by the dashed 
line in Fig. [5J The tail of the photodissociation cross sec- 
tion matches well with the one in the sequential picture 
(solid line) at about 0.2 MeV, as seen in the inner part of 
the figure. 

As we can see, for three-body energies smaller than 
0.2 MeV the dashed and solid curves are very different. 
Although still small in absolute values, the cross section 
in the direct picture (dashed line) can be orders of mag- 
nitude bigger than the one in the sequential picture (solid 
line). This difference is necessarily giving rise to different 
reaction rates at small temperatures. How sizable the dif- 
ference is between the reaction rates is the topic of the 
following subsection. 



4.2 Reaction rates 

In Fig. [7] the solid and open circles are the reaction rates 
for the triple alpha reaction given by Fowler et al. in 
Ref.[15j and the NACRE compilation [TO], respectively. 
The main difference between these rates is found at high 
temperatures, where the result from |15] is below the one 
of NACRE HBJ. This is due to the fact that while the ef- 
fects of the first 2 + resonance are included in [16], they are 
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Fig. 7. (color online) Reaction rate for the triple a process. 
The solid and open circles are the results from [15] (Fowler 
1967) and [16] (Angulo 1999), respectively. The thick solid 
curve is the calculation as in Eq.Q assuming a direct cap- 
ture process at very low energies. The thin solid line is the 
calculation in the sequential picture {TO]). The dashed line has 
been obtained in the extreme sequential approximation (|12[) . 
The inset shows the ratio between the thick solid rate (direct 
capture assumption at low energies) and the thin solid line 
(sequential capture assumption at low energies). 
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Table 1. Computed reaction rates for the a + a + a — > 12 C + 7 
reaction assuming a direct capture for the low energy tail of 
<t 7 (second column) and assuming a fully sequential process 
(third column). The fourth column gives the adopted reaction 
rate given in the NACRE compilation 16 . 



omitted in [15] . At low temperatures, both calculations, 
assuming both a fully sequential capture mechanism, pro- 
vide essentially the same result. In the calculations shown 
in this section the lowest 2 + resonance in 12 C has been 
included using the resonance parameters given in [16] . 

In the figure, the thin dashed line gives the reaction 
rate obtained in the extreme sequential model, i.e., when 
Eq.Q (solid curve in Fig.© is inserted into Eq.([T2]). This 
approximation amounts to cutting the tail of the cross 
section for energies smaller than the one of the interme- 
diate + state in 8 Be. This strong approximation gives 
rise to a reaction rate in good agreement with the results 
in |15] and |16] at high temperatures, while for temper- 
atures smaller than ~ 0.03 GK it clearly underestimates 
the rate. This is consistent with Fig. [3] where we see that 
a decreasing width of the three-body resonance increases 
the range of temperatures where the sequential and the 
extreme sequential calculations agree. 

A full sequential description of the process requires the 
use of Ea. lfTO"]) . Compared to the extreme sequential limit, 
this calculation includes now the low energy tails of the 
cross sections (j3]) and ©, which are given by (J5JJ) and 
([20]) . respectively. This is done in [16], and since we are 
using the same resonance parameters, we of course recover 
the same result (thin solid line in the figure). 

However, when the low energy tail in a 1 is assumed to 
be given by (|2"l)l . which corresponds to a direct capture 
(dashed curve in Fig. [5]), and the a 1 cross section ^ is 
inserted into Eq. , we then get the reaction rate shown 
by the thick solid line in Fig. [7] As we can see, for temper- 
atures smaller than ~ 0.07 GK, the reaction rate obtained 
assuming a direct capture at very low energies is several 



orders of magnitude bigger than when the sequential cap- 
ture is assumed. This is appreciated more quantitatively 
in the inset of the figure, where we show the ratio be- 
tween both reaction rates, direct and sequential. We can 
immediately see that the ratio increases when decreasing 
the temperature. For a temperature of 10 -3 GK we have 
obtained a reaction rate for the direct capture about 12 
orders of magnitude bigger than in the sequential picture. 
This difference reduces to about 7 orders of magnitude for 
T=0.01 GK. 

In table [JJ we give the computed reaction rates in the 
direct picture (second column) and in the sequential pic- 
ture (third column) for temperatures from 10 -3 GK to 
0.1 GK. In Ref. [16] a lower limit, a higher limit, and an 
adopted value are given for this reaction rate. Our results 
in the sequential approach lie in between the two limits 
given in 16 , although closer to the lower one. In the last 
column of table [TJ we give the reaction rates quoted in 
Ref. [16] as adopted rate. 

It is important to note that the computed rates are 
the ones obtained in the limiting cases of a fully direct 
or a fully sequential description in the very low energy 
region. If both processes compete, the computed reaction 
rate would then be found in between the thin and thick 
solid curves in Fig. which can be taken as the upper and 
lower limits to the true reaction rate. In any case, even if 
the process is considered to be fully direct, the increase 
in the reaction rate compared to the NACRE result at 
T=0.01 GK is of only 7 orders of magnitude (see table [TJ), 
and thus far smaller then the 20 orders of magnitude ob- 
tained in [34] . 
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We emphasize that the calculation in the direct picture 
has been made with the same value of E 1 (= 3.7 x 10 -3 
eV) as in the sequential picture. A change in _T 7 implies 
precisely the same change in the reaction rate, as seen 
immediately from Eqs.((3]) and ([3]). 



5 The 9 Be rate 

Let us consider the a + a + n — > 9 Bc + 7 reaction where the 
principles are similar to the case of 12 C. The properties of 
the low energy spectrum of 9 Be will determine the reaction 
rate at low temperatures. As for 12 C, the presence of the 
internal two-body subsystem, 8 Be with the very narrow + 
resonance at 0.092 MeV, suggests a sequential description 
of the capture process as the most appropriate model [161 

The main difference compared to 12 C is that now 9 Be 
does not show a low-lying and narrow three-body reso- 
nance like the Hoyle state. Instead, 9 Be has a low ly- 
ing l/2 + resonance in the vicinity of 0.11 MeV [37], only 
slightly above the 0.092 MeV of the + resonance in 8 Be. 
The width of the 1/2+ state is estimated to be around 
0.1 MeV, although 0.2 MeV is not numerically excluded 
[55] . As a consequence the photodissociation cross section 
shows a relatively broad peak at a three-body energy of 
around 0.11 MeV, in such a way that tr 7 is not negligible 
in the vicinity of the two-body resonance energy i? r =0.092 
MeV. The system under investigation is now similar to the 
schematic case shown in the lower part of Fig. SJ 



5.1 The photodissociation cross section 

The experimental photodissociation cross section of 9 Be 
can be found in [29], and it is shown in Fig. [8] by the solid 
circles as a function of the three-body energy E. Very re- 
cently new data for the peak corresponding to the 1/2+ 
resonance in 9 Be have been published in [39 , and they are 
shown by the open circles in the figure. The experimental 
data show that for three-body energies below the + res- 
onance at 0.092 MeV in 8 Be, the cross section essentially 
vanishes, or at least it is extremely small. The energy of 
the two-body + resonance is indicated by the arrow on 
the x-axis. 

This fact supports the assumption that the interme- 
diate + state in 8 Be is actually populated in the pro- 
cess, and therefore the sequential description appears to 
be appropriate. In fact, the energy dependence of the cross 
section © for sequential decay of 9 Be into 8 Be plus a neu- 
tron is governed by (|18p . which implies that <r 7 vanishes 
when E' = 0, i.e., when E = E r . Furthermore, for high 
temperatures, a pure three-body calculation of the reac- 
tion rate, without any additional assumption about the 
reaction mechanism, agrees reasonably well with the one 
obtained in the sequential picture [12) . 

As for 12 C, the behavior of ct 7 at very low energies, 
where again ct 7 takes very small values, will determine the 
value of the reaction rate at very low temperatures, and 
this rate could change depending on which model, direct or 
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Fig. 8. (color online) Photodissociation cross section for 9 Be 
as a function of the three-body energy E. The experimental 
data are from [29] (Sumiyoshi 2002, solid circles) and from 
[39] (Burda 2010, open circles). The dashed, dot-dashed, and 
dotted curves correspond to the contribution of the l/2 + ,5/2 + , 
and l/2~ states in 9 Be, respectively. The thick solid line gives 
the sum of the three contributions. The 5/2 + resonance is fitted 
as in [29] . The solid thin line is the total cross section when the 
5/2 + resonance is parametrized as in [16] . The vertical arrow 
in the x-axis indicates the energy of the + resonance in 8 Be. 
The inner part is a zoom of the low energy region of the cross 
section, where the solid and dashed lines are the cross section 
when the sequential and direct pictures are assumed in the low 
energy region, respectively. 



sequential, is assumed for the capture mechanism at such 
low energies. To investigate this issue we proceed as for 
12 C. At high temperatures the three-body and sequential 
pictures provide similar results (as we can see in the lower 
part of Fig. 2]). For energies higher than E r we then start 
by taking the simple parametrization of the cross section 
used in the sequential description. In particular, we again 
choose the resonance parameters given in [16] . and the 
1/2+ , 5/2 + , l/2~ resonances in 9 Be are included in the 
calculation. This gives rise to the cross section shown in 
Fig. |8j by the thin solid line, where the energy dependence 
given by Eq. fT?]) with A = 1, and Eq. ([T3j) . has been used 
for all the resonances. 

The peak corresponding to the 5/2 + resonance is clearly 
overestimated. To solve this problem, we have taken for 
this particular 9 Be resonance the energy and width used 
in [29] , which gives a much better agreement with the ex- 
perimental cross section, as shown by the thick solid line. 
In any case, the properties of the 5/2 + resonance play a 
minor role in the low temperature behavior of the reaction 
rate. The relevant resonance in this temperature region is 
the l/2 + state, for which similar parameters are employed 
in [TB] and [55]. In Fig.[3Jthe dashed, dot-dashed, and dot- 
ted curves show the contribution to the total cross section 
(thick solid curve) from the l/2 + , 5/2 + , and l/2~ reso- 
nances, respectively. 

In the inner part of Fig. [5J we show a zoom of the 
low energy region of the photodissociation cross section. 
When the sequential picture is assumed, the cross section 
takes the form ^ with E a b, c given by (|18l) . and then er 7 
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Fig. 9. (color online) Reaction rate for the reaction a+a+n — > 
9 Be + 7. The solid and open circles are the results from [29] 
(Sumiyoshi 2002) and [IB] (Angulo 1999), respectively. The 
dashed line has been obtained in the extreme sequential ap- 
proximation (|12|) . The thin solid line is the calculation in the 
sequential picture (JTOj) . The thick solid curve is the calcula- 
tion Q assuming a direct capture process at very low energies 
(dashed curve in the inset of Fig. [5]). The inset shows the ratio 
between the thick solid rate (direct capture assumption at low 
energies) and the thin solid line (sequential capture assumption 
at low energies) . 
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Table 2. Computed reaction rates for the a + a + n — > 9 Be + 7 
reaction assuming a direct capture for the low energy tail of 
<t 7 (second column) and assuming a fully sequential process 
(third column). The fourth column gives the adopted reaction 
rate given in the NACRE compilation 16 . 



vanishes for E = E r , as shown by the solid line in the inset 
of the figure. However, if for such low energies we assume 
a direct mechanism, due to the Coulomb repulsion of the 
two alpha particles, the low energy behavior of the cross 
section takes the form and the cross section behaves 
as shown by the dashed line in the inset of Fig. [8] 



5.2 Reaction rates 

The solid and open circles in Fig. [9] are the reaction rates 
for the a + a + n — > 9 Bc + 7 process given by Sumiyoshi 
et al. in Ref. [55] and the NACRE compilation [T£], re- 
spectively. They basically agree in the whole temperature 
range shown in the figure, although for low temperatures 
the results from [53] are slightly below the ones in [IB] . 
This is due to the fact that the treatment of the low energy 
tails in the cross section is not exactly the same in both 
cases. More precisely, for the 5/2 + resonance the widths 
involved in Eqs.© and (0 are taken as constant in [TB] . 
while in [29] they are assumed to be energy dependent. 
In the present work we use the same energy dependence 
as for the l/2 + resonance, and therefore our computed 
reaction rates in the sequential picture at very low tem- 
peratures should also be a bit smaller than the ones given 
in [IB]. 

If we now take the cross section given by the thick 
solid line in Fig. [S] and use it in Eq. (IT21 , we then get the 
reaction rate in the extreme sequential picture, which is 
shown by the dashed line in Fig. [9] As expected, this crude 
approximation fails completely at low temperatures. As we 
already know, a correct sequential calculation requires use 



of Eq. (fTU]) , which implies use of the cross section ^ and 
consequently inclusion of the energy tail ([2"T]) coming from 
the decay of the + resonance in 8 Be. When this is done 
we obtain the rate given by the thin solid line in Fig. [S] 
which, as expected, are slightly below the results given in 
[16] (open circles) at very low temperatures. 

Finally, let us assume that the low energy tail in 07 
corresponds to a direct decay mechanism (dashed curve in 
the inset of Fig. [5]). In this case, the cross section is given 
by ([!]), and the low energy tail takes again the form (J23|. 
The reaction rate is then given by Eq.®, and we obtain 
the thick solid curve in Fig. [3] The reaction rate computed 
assuming a direct capture at very low energies and the 
one obtained assuming a sequential capture begin to differ 
for temperatures smaller than about 0.07 GK. Very soon 
the rate obtained in the direct picture is several orders of 
magnitude bigger than the sequential one. This is better 
appreciated in the inset of the figure, where we show the 
ratio between both reaction rates. For a temperature of 
0.01 GK the reaction rate in the direct picture is almost 7 
orders of magnitude bigger than in the sequential picture. 
The general behavior of the reaction rates is similar to the 
one obtained for the triple alpha reaction, although due 
to the smaller Coulomb repulsion, the absolute values of 
the rates are now much bigger. 

As for 12 C, we now give in tabled] the computed re- 
actions rates in the direct picture (second column) and 
in the sequential picture (third column) for temperatures 
from 10~ 3 GK to 0.1 GK. In the fourth column we give 
the adopted value quoted in [16]. As already mentioned, 
for very low temperatures, our computed rates in the se- 
quential picture are below (up to more than a factor of 
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4 for T = 0.001 GK) the rates given in [16 . In fact, our 
values are even below the lower limit for the rate given in 
[IT)] . When a constant value is taken for the different -T's 
involved in the decay of the 5/2 + resonance, a rate within 
the limits given in |16j is obtained. 

The comments made when discussing the reaction rates 
for the triple alpha reaction are still valid here. First, the 
computed rates are the limiting cases of a fully direct pro- 
cess and a fully sequential process in the very low energy 
region. If both processes compete, the computed reaction 
rate would be found in between the thin and thick solid 
curves in Fig. O And second, the calculation in the di- 
rect picture has been made assuming that r 7 (E), given 
by Eq. (fTT|) . is the same as in the sequential picture, with 
jT 7 = 0.51 eV. A change in 7^ 7 would imply the same 
change in the reaction rate. 

6 Summary and conclusions 

In this work we investigate how different descriptions of 
the photodissociation cross sections at very low energies 
change the reaction rates for radiative three-body cap- 
ture processes. We focus on the reaction rates at the very 
low temperatures relevant for the nucleosynthesis of ele- 
ments in the core of a star. More precisely, we consider 
the radiative capture of three particles (nuclei, neutrons, 
protons. . . ) into a bound nucleus plus a photon. We inves- 
tigate how the reaction rate for such a process changes at 
low temperatures depending on the capture mechanism, 
i.e. either sequentially through an intermediate two-body 
state or directly. 

We first establish notation and definitions, and in par- 
ticular we specify the formulae for sequential and direct 
reaction rates for three arbitrary particles. These expres- 
sions rely heavily on the energy dependence of the pho- 
todissociation cross section. We have previously computed 
rates and cross sections in a full three-body computation 
without assumptions of reaction mechanism. In these com- 
putations the shape of the cross sections around resonance 
positions are computed as well, except at very low energies 
where the discretization is insufficient to characterize the 
resonance shapes, and a specific form has to be chosen. 
These low energies are crucial in the present work where 
we assume a simple Breit-Wigner form for the photodisso- 
ciation cross section but with energy dependent width as 
in i?-matrix analyses. This is not a severe limitation since 
the overall conclusions are general and very robust. 

We compute rates as functions of the involved pa- 
rameters for both two- and three-body resonances. We 
study the dependence on reaction mechanisms for differ- 
ent resonance positions and widths for different combina- 
tions of charges. The different rates almost coincide for 
temperatures somewhat larger than the three-body reso- 
nance energy. The direct decay only depends on position 
and width of the three-body resonance entering through 
the photodissociation cross section, and thus not on any 
two-body substructure. On the other hand, the sequential 
decay mechanism strongly depends on the relative posi- 
tion of the intermediate two-body resonance. If the two- 



body resonance is above the three-body resonance and the 
widths are substantially smaller than the energy difference 
between these resonances the sequential rate reproduces 
almost the extreme limit of zero widths where the rate 
expressions are very simple. 

When the two and three-body energies exchange posi- 
tions much larger rates are found and the sequential and 
direct rates apparently becomes more similar. However, 
they still deviate by many orders of magnitude strongly 
depending on charge and, without Coulomb potentials, 
angular momentum. In all cases the sequential rate is sub- 
stantially smaller than the direct rate, except when the 
two and three-body energies are fairly close to each other 
and all particles are charged. 

We turn to practical and realistic estimates of astro- 
physical reaction rates. We compute rates for the a + a + 
a — > 12 C + 7 and the a + a + n — > 9 Be+7 reactions. Three- 
body calculations [12, show that at high temperatures the 
usual sequential description [TMll)ll2"5] is appropriate. We 
focus on the reaction rates at low temperatures, where due 
to the extremely small energies, clearly below the energy of 
the available intermediate two-body state, the sequential 
picture is questionable. Therefore, for high temperatures 
we take one of the sequential parametrizations of the pho- 
todissociation cross section available in the literature. For 
the low energy tail we consider the possibility of a direct 
process for which the penetration factors are known for 
both Coulomb and centrifugal barrier potentials. 

We find that a direct description of the low energy 
tail of the photodissociation cross section enhances the 
reaction rates at low temperatures by several orders of 
magnitude compared to the sequential description. For a 
temperature of 0.01 GK the difference is of about 7 orders 
of magnitude for the two reactions investigated. These re- 
action rates in the direct and sequential pictures are the 
upper and lower limits of the rates. Therefore the result 
obtained for the triple alpha reaction in |34| is not consis- 
tent with our calculation. 

In conclusion, comparison between sequential and di- 
rect capture mechanisms reveal large differences at tem- 
peratures below the three-body resonance energy. These 
temperatures are important for the capture processes pro- 
ducing light nuclei in aging stars. The traditional use of 
the sequential models is insufficient and the direct reaction 
mechanism should be employed for these temperatures. 
Usually this leads to a substantial enhancement of the 
rates, but the precise increase depends on the character- 
istics of the most important resonances. Thus, individual 
calculations are necessary for each process. 
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